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EVIDENCE FOR AN OVERLUMINOSITY OF THE VARIABLE STAR RR LYR, 
AND A REVISED DISTANCE TO THE LMC 

M. Catelan 1 and C. Cortes 1,2 

ApJ (Letters), in press 

ABSTRACT 

We use theoretical models to establish a tight relationship for the absolute magnitudes of RR Lyrae stars 
as a function of their periods and Stromgren pseudo-color cq = (u— v)q — (v— b)o. Applying this to RR Lyr, 
and comparing the result with the predicted average absolute magnitude for stars of similar metallicity from 
the same models, yields an overluminosity of 0.064 ±0.013 mag in Stromgren y (and thus similarly in V) 
for RR Lyr. Based on a revised value for RR Lyr's trigonometric parallax, and on a newly derived reddening 
value of E(B-V) = 0.015 ±0.020, we provide a corrected relationship between average absolute magnitude 
and metallicity for RR Lyrae stars that takes RR Lyr's evolutionary status fully into account for the first time. 
Applying this relationship to the LMC, we derive a revised true distance modulus of (m-M)o = 18.44 ±0.1 1. 

Subject headings: stars: distances — stars: horizontal-branch — stars: variables: other — distance scale 
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1. INTRODUCTION 

RR Lyrae (RRL) stars are radially pulsating variable stars 
with periods in the 0.2-1.0 d range, and visual ampli- 
tudes < 2 mag. They are accordingly easily identified, and 
play a key role as the cornerstone of the Population II dis- 
tance scale. They are extensively used to determine dis- 
tances to old and sufficiently metal-poor systems, where they 
are commonly found in large numbers. In particular, RRL 
are present in globular clusters (GC's) and the d warf galax- 
ies in the neighborhood of the Milky Way (e.g. jGreco et alj 
20071). and h ave also been identified in the M31 field (e.g., 



Brow n et alll2004t iDolphin et alj |2004). in some M31 com- 
panions (e.g.JPritzl et alj|2005l) . and in at least 4 M31 GC's 



dClementini et al.ll2001l) T 

While period-luminosity (PL) relations for 
stars in JHK ha ve been known and studied 
long time (e.g.. lLongmore. F ernlev. & Jameson 
Cate lan. Pritzl. & Smithl 12004 Del Principe etalj 



RRL 

for a 



119861 
l2006h . 

only recently have s uch relations bee n exten ded to bluer 
filters. In particular, ICortes & Catelan! d2007l) have carried 
out a theoretical study of RRL stars in the lStromgrenl (119631) 
uvby passband system, concluding that several relationships 
involving the (rather reddening-independent) pseudo-color 
Co = («— v)o — (v— b)o are present that may prove helpful 
in deriving distances to RRL stars. As well known, the 
Stromgren system represents an invaluable tool in the study 
of the physical parameters of stars, including tempera- 
ture, surface gravity, metallicity, and (clus ter) ages (e.g., 
ISchuster & Nissedll989tlGrandahl et alj|2000h . 

In this Letter we provide, for the first time, a method to 
precisely estimate the degree of overluminosity in the visual 
for even field RRL stars. This is based on a comparison be- 
tween the predicted average RRL absolute magnitude for the 
star's metallicity, on the one hand, and its own precise abso- 
lute magnitude, on the other - the latter being derived based 
solely on cq and the pulsation period. 

2. THE METHOD 
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The chemical c ompo sition of RR Lyr was studied by 
IClementini etail (119951) . who derived [ Fe/H| = -1.39 
[a/_Fe] = +0.31 for the star. Accordi ng to iBragaglia et all 
feoOll) . this [Fe/H] value is in the IZinn & West! (11984 
ZW84) scale; thi s transl ates to [Fe/H] = -1.16 in the 
ICarretta & Grattonl (1 19971 CG 97) scale. Based on 



71 

iSalaris. Chieffi. & Stranierol (Q993), we estimate a metallic- 
ity Z = 1.4 x 10~ 3 or 2.3 x 10~ 3 in the two scales, respec- 
tively. Our results will thus be based on models computed 
for a metallicity Z = 2 x 10~ 3 . 

We have computed extensive sets of horizontal branch (HB) 
simulations following the same procedures already des cribed 
elsewhere dCatelan et al.ll2004tlCortes & Catelanll2007l) . The 
main sequence helium abundance in our models is 1ms = 0.23. 
In Figure [TJ we show a synthetic M y , b-y CMD, for a rather 
even HB type. The blue e dge of the instabilit y strip (IS) 
is based on equation (1) in Capu to et al.l d 1987b . shifted in 
temperature by -200 K. The IS width was been taken as 
Alog T e ff = 0.075; this defines the IS red edge. These choices 
provide good agreement with recent theoretic al prescriptions 
and the observations (see §6 in ICateladl2004l) . 

RRL-based distances are usually derived using a simple re- 
lation between average RRL V-band magnitude and metal- 
licity. However, RR Lyr is the only RRL-type star with a 
sufficiently a ccurate trigonometric parallax for distance de- 
terminations dBenedict et al.l [2002). This star has accord- 
ingly been used by several authors to place the zero point 
of the (My) - [Fe/Hl relat ion on a firmer footing (e.g., 
ICacciari & Clem entini 2003). Unfortunately, RR Lyr's evo- 
lutionary status remains at present unknown, and therefore it 
is still unclear whether it is representative of the underlying 
stellar population to which it is associated. This is an impor- 
tant issue since, as can be seen from Figure [T] a scatter in 
RR L magnitudes (by up to a few tenths of a magnitude; see 
also[Sandage 1990) is always expected within any single pop- 
ulation, and therefore distances based on current calibrations 
of th e (My)- [Fe/H l relation may be correspondingly in error 



(e.g.. lCatelan1l2005l) . 

To emphasize this point, Figure [2] compares the difference 
in absolute magnitude between individual RRL stars in three 
HB simulations with widely different HB morphologies and 
the average HB absolute magnitude computed for a large sam- 
ple of RRL stars in extensive HB simulations covering the full 
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FIG. 1. — Simulated CMD for an intermediate HB type and a chemical 
composition Fms = 0.23, Z = 0.002. Non-variable HB stars are shown as 
light gray circles, whereas RRL variables are indicated as dark gray circles. 
The position of RR Lyr is shown as a red circle with error bars (see text). 
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FIG. 3. — As in Figure[2] but comparing the absolute magnitudes provided 
by equation (TJ, here shown as M v (fit), and the absolute magnitudes of each 
individual star M,,(sim), as provided by the HB simulations. Note that the 
scatter is at the level of a = 7.2 X 10~ 3 mag only. 
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FIG. 2. — Difference between each RRL star's individual My value and the 
average over all 8114 synthetic RRL stars in our simulations, for a blue HB 
simulation (blue diamond signs), an even HB simulation (gray circles) and a 
red HB simulation (red squares). 

range of observed HB morphologies ((M y ) = 0.6614 ±0.0008 
and median M y = 0.6896, based on 81 14 synthetic stars). 

Up to now, no method had been available to infer the evo- 
lutionary status of individual RRL stars to a precision better 
than 0.1 mag. The Stromgren system may provide us with 
a solu tion to this problem. In particular, Cortes & Catelanl 
(2007) have shown that, by incorporating a pseudo-color term 
to the PL relation, exceedingly tight relations derive, even 
for the bluer bandpasses. We have thus used the sample of 
8114 synthetic RRL stars mentioned above to search for high- 
precision analytical fits that would provide absolute magni- 
tudes to a precision of ~ 0.01 mag. We thus found a linear 
dependence on logP and a cubic dependence on the (log of 
the) pseudo-color c . Our final relations thus read as follows: 



and 



: a + b In c + c (In c ) 2 + d (In c ) 3 + e log P. 



(D 



The coefficients of the fits for each of the four uvby bands, 
along with their respective errors, are given in Table [T] The 
correlation coefficient ranges from r = 0.992 to 0.998, and the 
standard errors of the estimates from 0.0072 to 0.0078 mag. 

To check how well equation (HJ is able to reproduce M y 
for individual RRL stars, we show, in Figure [3] the differ- 
ence between the M y value implied by this equation and the 
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Coefficients of the Fits 
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+1.8000 


0.0010 
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+0.4535 


0.0009 
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-1.1484 


0.0029 
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-1.9103 


0.0027 
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-0.4903 


0.0213 
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-0.9250 


0.0204 
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-1.3840 


0.0463 
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0.0443 
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0.0026 
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-1.8575 


0.0025 
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+0.1070 


0.0009 
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-0.1244 


0.0009 
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-1.6721 


0.0027 
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-1.3546 


0.0026 
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-0.8628 


0.0197 
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-0.7413 


0.0195 
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-1.2972 


0.0430 
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-1.3146 


0.0425 
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-1.9507 


0.0024 


e 


-2.0577 


0.0024 



input value from the same HB simulations as in Figure [2] As 
can be clearly seen, the scatter, which in the previous figure 
amounted to several tenths of a magnitude, has been reduced 
to a few thousandths of a magnitude only. It is thus clear that 
equation (Q~|i is capable of providing M v values for RRL stars 
with remarkable precision, provided of course their periods 
and Co values are known. Comparing this with (M y ) from the 
same models, one immediately obtains a robust estimate of 
the degree of overluminosity for individual RRL stars. 

3. A NEW LOOK AT RR LYR 

Though RR Lyr is an extensively studied star, not many 
RRL surveys have bee n carried out in the Stromgren sys tem 
to date - the papers bv lEpsteinl (fl969l) and lSieie1(fT982l) be- 
ing notable exceptions. We shall adopt the latter' s dataset in 
this work. ISiegell carried out his uvbyf3 study of RR Lyr using 
the 16 in telescope at Kitt Peak, and provided a c\ light curve. 
This is shown in FigureHl where Siegel's two datasets (repre- 
senting two different phases in the star's ~ 40 d-long Blazhko 
cycle) are displayed. 

Naturally, in order to be able to apply our method to RR 
Lyr, we first need to compute a suitable average c\ value for 
the star. In this sense, we have first carried out a Fourier fit to 
the star's observed c\ curve. The result of a 5th-order Fourier 
decomposition is also shown in Figure|4] Averaging over this 
curve in magnitude units, we find (ci) mag = 0.861 ±0.002. 
Performing the same average in intensity units, we obtain a 
slightly different result, namely (ci) = 0.850 ± 0.002. In like 
vein, (u) -2(v) + (b) = 0.850 ± 0.010. We have performed the 
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FIG. 4. — The ci data for RR Lyr (from ISieeell 19821) are plotted as a 
function of the pulsation phase. The two datasets from Siegel are given with 
different symbols. A 5th-order Fourier fit to the data is also shown. 

same test using Fourier decompositions of different orders, 
with no significant differences in the final average values. 
Likewise, the average c\ value, as derived based on Siegel's 
(1982) dataset 2 only, is not noticeably different from those 
provided above. (His dataset 1 is missing the minimum.) One 
of the main advantages of the present method is indeed the 
fact that the amplitude of the c\ light curve is rather small 
(0.52 mag, according to the 5th-order Fourier fit), so that the 
differences between the several procedures to obtain the av- 
erage quantities appear to be minor. We finally adopt for RR 
Lyr an average c\ = 0.855 ± 0.005. 

As far as the reddening correction goes, we remind the 
reader that c\ is a rather reddening-insensitive index, with 
co = c\ -0.1 5E(B-V). The redd ening of RR Lyr is not well 
determined. Sollima et al.| (120081) have recently argued that a 
value E(B-V) = 0.02 ±0.03 safely covers all of the available 
determinations. An independent reddening determination is 
afforded by comparing the star's Stromgren photometry and 
our relations (see eq.Q]and Table 1). 

In this sense, we have computed magnitude- and intensity- 
averaged b—y, v-b, u — v, and u—y colors for RR Lyr, and 
compared the values derived using Siegel's (1982) combined 
datasets and his dataset 2 only. We then averaged the results, 
adopting as final average colors for RR Lyr the following val- 
ues: b-y = 0.282 ±0.002, u-v= 1 .215 ± 0.003, v-b = 0.357 ± 
0.002, and u-y = 1.857 ±0.009. The error bars include the 
errors of the Fourier fits used to compute the means, as well 
as the dispersion around the mean value given by the differ- 
ent definitions of the "equivalent static star." Comparing these 
values with the intrinsic colors implied by our theoretical re- 
lations at c = 855 ± 0.005 and a period P = 0.5668386 d 
dKolenberg et all 120061) . using the relationships E(b—y) = 
0.74E(B-V), E(u-v) = 0.65 E(B-V), E(v-b) = 0. 50 E(B-V) 
and E(u-y) = U >9E(B-V) dCrawford & Mandwewd jfl97a 
IClem et alj2004l) . and taking a weighted average over the thus 
derived individual E(B-V) values, we find for RR Lyr a re- 
vised reddening value of 

E(B-V) = 0.015 ±0.020, (2) 

which also implies a small correction to the adopted cq value 
for RR Lyr, which becomes c = 0.853 ± 0.006. 

Inserting these values into equation ([TJ for y, we find M y = 
0.585 ±0.010 mag. Using the average M y for RRL stars of 
similar metallicity, from §2, we find that RR Lyr is an overlu- 



minous star, by about 0.077 ±0.010 mag in y (and thus V)? 
The inferred position of the star is shown in the CMD of Fig- 
ure Q] For an overluminosity in y of zero to obtain, we would 
have needed a cq = 0.801 for RR Lyr, which is many sigma 
away from the value derived on the basis of Siegel's (1982) 
data. Note that the star's overluminosity, with respect to the 
median, is even higher: 0.105 ±0.010 mag. We emphasize 
that this result is based on a differential use of the theoretical 
models. Within the error bars, RR Lyr's overlumi nosity is also 
not inc onsistent with the recent near-IR study bv lSollima et al.l 
(2008, their Fig. 2). 

To check the model d ependence of this result, we repeated 
our procedur e using the iPietrinferni et ail (|2004) tracks, cor- 
rected as in ICassisi et al.l (120041) ! Based on the resulting 
HB simulations we rederive the parameters of equation ([T), 
leading to a relation with r = 0.993 and standard error of 
the estimate 0.006 mag. We infer the M y value for RR Lyr 
therefrom (M y = 0.580 mag), and compare with the value 
(My) = 0.6576 mag from the same simulations. This gives for 
RR Lyr an overluminosity of 0.077 mag, in p erfect agreement 
with t he value derived earlier. Using the iLee & Demarquel 
( 1990) models for Ims = 0.23, the overluminosity we find is 
0. 064 mag, again consiste nt with our previous results. 

ICorte s & Catelan (2007) caution that a star's exact overlu- 
minosity will depend on the adopted metallicity scale. To 
safely accomodate the resulting range of possible values for 
RR Lyr, we finally adopt for this star an overluminosity in y 
(and thus V) of 0.064 ± 0.013 mag. 

4. A REVISED RRL (Afy) - [FE/H] RELATION 

Let us assume, as often done, a relation of the form (My) = 
a[Fe/H] + (3. Here we shall adopt a slope a = 0.23 ± 0.04, 
as suggested in several rece nt reviews of the subject (e.g., 
Cacciari & Clementini 2003). What about the zero point (31 

iBenedict et alj d2002l) obtained, using the Hubble Space 
Telescope (HST), an accurate value for RR Lyr's trigonomet- 
ric parallax, namely 7r^ T = 3.82 ± 0.20 mas. Very recently, 
Ivan Leeuwenl (12007) revised the trigonometric parallaxes pro- 
vided by Hipparcos, and arrived at a value 7r ab ^ = 3.46 ± 
0.64 m as for RR Lyr. A weighted avera ge of ground-based 
studies dvan Altena. Lee. & Hoffleiilll995l) indicates a paral- 
lax 7rfb° Und = 3.0 ± 1.9 mas for this star. Taking a weighted 
average of these results, we obtain a final value of 7r a b s = 
3.78 ± 0.19 mas for RR Lyr. This implies a revised distance 
modulus of (m-M)o = 7.11 ±0.11 mag for the star. 

In recent wor k, an intensity-weig hted mean magnitude of 
(V) = 7.76 mag dFernlev et al.lll998l F98) has been adopted 
for RR Lyr. However, we note that this value is based on 
Hipparcos photometry, which may require a non-trivia l trans- 
forma tion to the standard system. For comparison. iLavdenl 
d 19941) determines a (V) = 7.66 m ag. andlLavden et all dl996h 
find instead (V) = 7.74 mag. iGould & Popowskil (U998l 
GP98) argue strongly in favor of the Hipparcos-based mag- 
nitudes of F98, but propose a (small) reddening-related cor- 
rection: Vgp98 = Vf98-0.2£ , (B-V). Using our derived redden- 
ing value (eq. [2]), and adopting a standard extinction law with 
A v = 3.1 E(B-V), one finds Aj* RL = 0.045 ±0.062 mag. The 
final, extinction-corrected RR Lyr mean magnitude is accord- 
ingly (V ) =7.710 ±0.062 mag. 

3 To derive an overluminosity in other bandpasses is more tricky, in view 
of the non-zero slope of the HB in planes other than My, b — y. Taking the 
HB slope into due account, we find for RR Lyr overluminosities in u, v, b of 
-0.075, -0.113, and -0.096 mag, respectively. 
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As well known, the intensity- or magnitude-weighted mean 
magnitude does not necessarily correspond to the magnitude 
of the "equivalent static star" (i.e., the magnitude the star 
would have if it were not pulsating): an amplitude-dependent 
correction has to be applied. Such a cor r ection has been ob- 
tained by iBono. Caputo. & StellingwerS (1 19951) on the basis 
of detailed hydrodynami cal models. Inspection of the light 
curve s for RR Lvr (e.g.. ISzeidl & Kolllthl l2000t ISmith et alj 
120031; Ikolenberg et al. 200g) suggests that the amplitude in V 
oscillates in the 0.5 -1.1 mag range. According to Table 2 in 
IBono et all this is precisely the amplitude range over which 
the intensity-weighted mean magnitude provides the most ac- 
curate description of the magnitude of the equivalente static 
star. Taking, accordingly, a value (Vb) = 7.710 ±0.062 mag 
and a distance modulus (m-M)o = 7.11 ±0.11 mag, one finds 
an absolute magnitude for the star of My = 0.600 ±0. 126 mag. 

This, as previously argued, should not be viewed as rep- 
resentative of the average absolute magnitude of RRL stars 
of similar metallicity: according to our results, the latter are 
fainter (on average) than RRL by about 0.064 ±0.013 mag. 
Therefore, the average absolute magnitude of RRL stars of 
metallicity similar to RR Lyr's should be 



(M y ) = 0.664 ±0.127 mag. 



(3) 



Adopting for RR Lyr the metallicity values provided in §2, 
one then finds the following final relations: 

{M v )= (0.23±0.04)[Fe/H] ZW 84+(0.984±0.127), (4a) 
(M v )= (0.23 ±0.04)[Fe/H] CG 97+(0.931 ±0.127). (4b) 

To close, we note, in passing, that the result given in equa- 
tion (f3]i is in reasonable agreement with the average M y value 
at Z = 0.002 derived from our HB simulations (§2). 



5. A REVISED LMC DISTANCE MODULUS 

Equation (|4a]l implies an My = 0.644 ±0.141 mag at 
[Fe/H] = -1.48 ±0.07, which is the mean metalli city de- 
rived for LMC RRL variables by lGratton etaf] (120041 ZW84 
scale). Usin g a value (Vb) = 19.068 ±0.102 mag from 
iGratton et all one then finds a true distance modulus for the 
LMC of {m-M)\ wc = 18.42 ±0.17. If one uses instead the 
av erage values for LMC R RL stars independently determined 
bv lBorissova et all (|2004|) . namely [Fe /H] = -1 .46 ± 0.09 dex 
and (V) = 19.45 ± 0.04 mag, with their favored reddening of 
E(B-V) = 0.11 mag, one finds (V ) = 19.11 ±0.04 mag for 
a My =0.648 ±0.142 mag; thus {m-M)\ MC = 18.46 ±0.15. 
Taking a weighted average over these two results, we arrive at 
the following distance modulus for the LMC: 



{m-M)\ MC : 



18.44±0.11. 



(5) 



6. CONCLUSIONS 



In this Letter, we have established, based on Stromgren 
photometry, that the star RR Lyr is approximately 0.064 mag 
brighter than the average for other RRL stars of similar metal- 
licity. We also derive a new reddening value for RR Lyr, 
namely, E(B-V) = 0.015 ±0.020. These results, along with 
a newly derived trigonometric parallax for RR Lyr, lead to a 
revision in the RRL absolute magnitude-metallicity relation, 
and thus to a revised distance modulus for the LMC. 
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